Estimation of photon counting statistics with imperfect detectors
Han Xiao-Chuan1, Zhuang Dong-Wei1, Li Yu-Xuan1, Song Jun-Feng1, †, Zhang Yong-Sheng2, ‡
State Key Laboratory on Integrated Opto-electronics, College of Electronic Science and Engineering, Jilin University, Changchun 130012, China
Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, China

 

† Corresponding author. E-mail: songjf@jlu.edu.cn yshzhang@ustc.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 61627820, 11674306, 61590932, and 61377048).

Abstract

The study on photon counting statistics is of fundamental importance in quantum optics. We theoretically analyzed the imperfect detection of an arbitrary quantum state. We derived photon counting formulae for six typical quantum states (i.e., Fock, coherent, squeeze-vacuum, thermal, odd and even coherent states) with finite quantum efficiencies and dark counts based on multiple on/off detector arrays. We applied the formulae to the simulation of multiphoton number detections and obtained both the simulated and ideal photon number distributions of each state. A comparison between the results by using the fidelity and relative entropy was carried out to evaluate the detection scheme and help select detectors for different quantum states.

PACS: 42.50.-p
1. Introduction

Quantum optics is recognized as a counter-intuitive theory because of its non-classical properties. It has been widely used in applications such as quantum information technology,[15] quantum teleportation,[68] and quantum computation.[912] Photon number detection[13,14] is considered as an essential branch of quantum technology in photonic systems. The unknown quantum state to be measured does not always consist of only one photon. For multiphoton detections, photon-number-resolving (PNR) detectors[1517] are regarded as the simplest apparatus as the photon number distribution can be provided directly. However, they are not widely used because of their low efficiency.[18] The on/off detectors made of avalanche photodiodes exhibit a high efficiency,[19,20] but the number of incoming photons cannot be distinguished. Sperling et al. proposed a robust method for photon detection by using an on/off detector array instead of the PNR detector.[21] The schematic is illustrated in Fig. 1, where an unknown quantum state is divided equally into N ports via a series of hypothetical unitary operations. At each port, there is an on/off detector for photon detection. The total number of clicks of the detectors can then be counted, which is approximately recognized as the photon number distribution. The unitary operations here are assumed to be ideal.

Fig. 1. (color online) Schematic of multiphoton detections based on an on/off detector array. An unknown quantum state is equally divided into N ports by a series of hypothetical unitary operations (U(N)). At the end of each port, an on/off detector is ready for the detection of a possible incoming photon. The photon number distribution is obtained by summing the clicks of all the on/off detectors.

As it is possible that more than one photon can enter the same port in this method, the results may not be acceptable when there are insufficient on/off detectors. However, when there are sufficient on/off detectors, the probability of this situation is lower, and the detection results can then approach those of the ideal case. Here, we proposed a method using the fidelity and relative entropy to characterize the performance of the detectors. The photon counting formulae could be deduced based on multiple on/off detector arrays considering both the quantum efficiency and dark counts. Using these formulae, we took six typical quantum states as examples to study the number of detectors that can be used to obtain an ideal result as well as the effects arising from the noise counts.

2. Photon counting formulae based on on-off detectors

In the method based on the on/off detector arrays, the photon-counting formula of the quantum state |ψ⟩ can be written as[22]

where is the density operator of the quantum state |ψ⟩, N is the number of on/off detectors, and
is the projection operator.[21] Here, η is the quantum efficiency, υ is the dark counts, and : : is the normal ordering prescription. By combining Eqs. (1) and (2), we can obtain the detection probability when k detectors click (see Appendix A)
where
Here, S(m+l,k) is the Stirling number of the second kind,[23] is the binomial coefficient, and is the permutation coefficient. As only the term Ql(ψ, η) contains the information of quantum states and Ml(N,k,υ) is only the normalization coefficients, we can write Ql(ψ,η) in the form of different quantum states depending on their unique expressions as follows.

For the Fock state |n⟩,

For the thermal state,

For the squeezed-vacuum state,

where Pl(z) is the Legendre function,[24] A = [cosh2 ξ − (1 − η)2 sinh2 ξ]−1/2, x = A(1 − η)sinhξ, and y = sinhξ (see Appendix B). For a coherent state |α⟩ = exp(αa +α *a)|0⟩, the photon counting probability pk can be written as
For a coherent superposition state |ψ⟩ = M(|α⟩ + exp(iφ)| − α⟩), where |α⟩ is a coherent state, the photon counting probability pk can be written as

According to the above formulae, we calculated the photon number distribution for each state based on 4, 16, 64, and 256 on/off detectors respectively, and the results are shown in Fig. 2, which verify what we will illustrate in the last section.

Fig. 2. (color online) Six quantum states with a mean photon number of 3 detected by on/off detector arrays. Photon number distributions of (a) coherent, (b) thermal, (c) Fock, (d) squeeze-vacuum, (e) even coherent, and (f) odd coherent states with N = 4, 16, 64, and 256, respectively.
3. System similarity and fidelity

To compare the simulated photon number distributions with the ideal ones, we introduced a computable parameter, fidelity,[25] which is expressed as

where is the probability of the ideal photon number distribution and pk is the probability we calculated above. Equation (9) can also be understood as the similarity of two classical probability distributions. Assuming all detectors are perfect, we first calculated the fidelities of six quantum states with different numbers of on/off detectors, and the results are presented in Fig. 3(a). It can be clearly observed that as the number of detectors increases, the fidelity tends to be unity. For squeeze states, their photon number has a higher probability of 1 or 2, and the fidelity is fluctuant when less detectors are used.

Fig. 3. Six quantum states with a mean photon number of 3 detected by on/off detector arrays. (a) Calculated fidelities of each state based on perfect on/off detector arrays. (b) and (c) Calculated fidelities of each state based on 16 on/off detectors with finite quantum efficiency and dark counts. (d) Calculated fidelities of Fock state based on imperfect on/off detector arrays for ν = 5 × 10−6 and η = 0.9.

Figure 3(b) shows the variation of the fidelity as a function of quantum efficiency for the six quantum states. We find that with the quantum efficiency increasing, the fidelities of all the quantum states increase. In particular, the fidelities of the coherent and thermal states, which are close to the ideal states, change slowly, while the fidelities of the other four states change sharply.

This phenomenon is mainly due to the original photon number distributions of these quantum states. States that have non-continuous distributions, such as Fock, squeeze-vacuum, and odd coherent states, have intervals in their photon number distributions. After being detected by the on/off detectors, higher photon number distributions will partly transfer to lower ones. Therefore, the probability of the original photon number states is lower, while that of the adjacent photon number states is higher, resulting in a sharp change in fidelity. On the contrary, for states that have a continuous distribution, both the higher and lower photon number distributions will partly transfer to a lower photon number distribution. As a result, each photon number state will change by a small amount. In other words, the corresponding fidelity will change slowly. The processes are illustrated in Fig. 4.

Fig. 4. Photon number distributions of both ideal and simulated ones for (a) and (b) calculated squeeze-vacuum state, (c) and (d) calculated coherent state.

Figure 3(c) shows the variation of the fidelity as a function of dark counts. We find that all fidelities decrease as the dark counts increase except for the thermal state. Since dark counts are generated from thermal noise, it can compensate the same loss of photon counts due to finite detectors.

Figure 3(d) shows the variation of the fidelities of the Fock states as a function of the mean photon number and the number of detectors for ν = 500 s−1 × 10 ns = 5 × 10−6 (dark count v is the product of a single photon counter’s response time and dark counts per second) and η = 0.9. These parameters can easily be obtained in experiments. Besides the information we obtained in Fig. 3(a), we also find that all fidelities decrease as the mean photon number increases. This phenomenon is mainly due to the higher probability for more than one photon entering in the same port.

4. System diversity and relative entropy

In thermodynamics,[26] entropy is often used to characterize the degree of disorder of a system. In quantum systems, we can also use relative entropy to describe photon number distributions, which can be written as

Information entropy[27] is known to indicate the degree of ordering of a system.

It is quite different from the fidelity in the sense that the more ordered a system is, the higher the information entropy will be. Therefore, the relative entropy[28] can describe the distance between two probability distributions in some degree, which, in other words, can also be interpreted as the system diversity. A lower relative entropy here indicates that the simulated photon number distributions are approaching to the ideal distribution, in contrast to the fidelities. Figure 5(a) plots the calculated relative entropy of the six quantum states as a function of the number of on/off detectors. Like the fidelity, the relative entropy can also be used to distinguish between the photon number distributions simulated from the on/off detector array and the ideal distribution. We calculated the relative entropy with the increasing noise. The results are presented in Figs. 5(b) and 5(c), which show similar phenomena to the fidelities we illustrated in Fig. 4. We also calculated the variation of the relative entropy of the Fock states as a function of the mean photon number and the number of detectors for ν = 5 × 10−6 and η = 0.9, as shown in Fig. 5(d). It shows the same information which we have obtained in Fig. 3(d). In addition, the relative entropy can be used not only in quantum information, but also in classical information. In conclusion, the relative entropy may be more useful in future photon counting statistics.

Fig. 5. Six quantum states with a mean photon number of 3 detected by on/off detector arrays. (a) Calculated relative entropy for the six states as a function of detector number. (b) and (c) Calculated relative entropy of each state based on 16 on/off detectors with finite quantum efficiency and dark counts. (d) Calculated relative entropy of the Fock state based on imperfect on/off detector arrays for ν = 5 × 10−6 and η = 0.9.
5. Conclusion

We have derived a series of universal photon counting formulae for six common quantum states based on multiple on/off detectors. The formulae can be applied to obtain ideal and simulated photon number distributions. We have simulated the photon number distributions and found that increasing the number of detectors can result in a higher performance in both perfect and imperfect detections. We then derived and calculated the fidelity of the simulated photon counting distributions, which can assist us in choosing a suitable number of detectors to obtain a near-perfect photon number distribution of each state. We have also derived and calculated the relative entropy of the obtained photon counting distributions, which exhibit an opposite tendency to that of the fidelity. We found that for states with continuous distributions, the detection performance is better than that of the states with non-continuous distributions. Therefore, we should choose near-perfect detectors for non-continuous state detections, while imperfect detectors are acceptable for continuous state detections.

Reference
[1] Kervalishvili P J 2015 IEEE Seventh International Conference on Intelligent Computing and Information Systems 12–14 Dec 2015 Cairo, Egypt 1 10.1109/IntelCIS.2015.7397187
[2] Spiller T P 2003 Materials Today 6 30
[3] Erhard M Fickler R Krenn M Zeilinger A 2018 Light: Science & Applications 7 17146
[4] Caspani L Xiong C Eggleton B J Bajoni D Liscidini M Galli M Morandotti R Moss D J 2017 Light: Science & Applications 6 e17100
[5] Liu G Q Pan X Y 2018 Chin Phys. 27 020304
[6] Okada M Kan T Fuwa M Takeda S Yoshikawa J I Loock P V Furusawa A 2017 Conference on Lasers and Electro-Optics Europe & European Quantum Electronics Conference 25–29 Jun. 2017 Munich Germany 2017 1 1
[7] Djordjevic I B 2016 IEEE Photonics Journal 8 1
[8] Bouwmeester D Pan J W Mattle K Eibl M Weinfurter H Zeilinger A 1997 Nature 390 575
[9] Liang Z F Ma S L Xue H T Ding F Liu J L Tang F L 2018 Chin. Phys. 27 016801
[10] Lu F F Wang C F 2016 Chin. Phys. Lett. 33 74202
[11] Ren X Z He S Cong H L Wang X W 2012 Acta Phys. Sin. 61 124207
[12] Wakakuwa E Soeda A Murao M 2017 IEEE Transactions on Information Theory 1 1 10.1109/TIT.2017.2709754
[13] Helstrom C W 1969 Journal of Statistical Physics 1 231
[14] Zhang J Itzler M A Zbinden H Pan J W 2015 Light: Science & Applications 4 e286
[15] Munro W Nemoto K Beausoleil R Spiller T 2005 Phys. Rev. 71 033819
[16] Rosenberg D Lita A E Miller A J Nam S W 2005 Phys. Rev. 71 061803
[17] Rohde P P 2005 Journal of Optics B: Quantum and Semiclassical Optics 7 82
[18] Rohde P P Webb J G Huntington E H Ralph T C 2007 New J. Phys. 9 233
[19] Eisaman M Fan J Migdall A Polyakov S V 2011 Review of Scientific Instruments 82 071101
[20] Cova S Ghioni M Lacaita A Samori C Zappa F 1996 Appl. Opt. 35 1956
[21] Sperling J Vogel W Agarwal G S 2012 Phys. Rev. 85 023820
[22] Hall B C 2013 Quantum Theory for Mathematicians Berlin SpringerVerlag 53
[23] Adamchik V 1997 Journal of Computational and Applied Mathematics 79 119
[24] Abramowitz M Stegun I A 1972 Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (9th printing) Denver Dream Books Company, LLC 331
[25] Cozzini M Ionicioiu R Zanard P 2007 Phys. Rev. 76 104420
[26] Callen H B 1998 Thermodynamics and an Introduction to Thermostatistics 27
[27] Tan Y Z Wu C F 2003 Journal of Natural Resources 1 017 (in Chinese)
[28] Kullback S Leibler R A 1951 The Annals of Mathematical Statistics 22 79